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ABSTRACT 

We present in this paper a new three-dimensional galaxy classification system 
designed to account for the diversity of galaxy properties in the nearby universe. 
To construct this system we statistically analyse a sample of > 22, 000 galaxies at 
v < 15,000 km s _1 [z < 0.05) with Spearman rank and principal component 
analyses. Fourteen major galaxy properties are considered, including: Hubble type, 
size, colour, surface brightness, magnitude, stellar mass, internal velocities, HI gas 
content, and an index that measures dynamical disturbances. We find, to a high degree, 
that most galaxy properties are correlated, with in particular Hubble type, colour, 
and stellar mass all strongly related. We argue that this tight 3-way correlation is 
a result of evolutionary processes that depend on galaxy mass, as we show that the 
relation between colour and mass is independent of Hubble type. Various principal 
component analyses reveal that most of the variation in nearby galaxy properties can 
be accounted for by eigenvectors dominated by: (i) the scale of a galaxy, such as its 
stellar mass, (ii) the spectral type, and (iii) the degree of dynamical disturbances. We 
suggest that these three properties: mass, star formation, and interactions/mergers 
are the major features that determine a galaxy's physical state, and should be used to 
classify galaxies. As shown in Conselice (2003), these properties are measurable within 
the CAS (concentration, asymmetry, dumpiness) structural system, thus providing an 
efficient mechanism for classifying galaxies in optical light within a physical meaningful 
framework. We furthermore discuss the fraction and number density of galaxies in the 
nearby universe as a function of Hubble type, for comparison with higher redshift 
populations. 
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1 INTRODUCTION 

Galaxies as distinct stellar, gaseous, and dark matter units 
have dozens of properties measurable by a variety of tech- 
niques from gamma-ray to radio wavelengths. Some of the 
most basic properties include: spectral energy distributions, 
spectral line properties, sizes, internal velocities, structural 
features, and morphologies. From these observable quanti- 
ties we can derive additional physical information, such as 
stellar, gaseous, and total masses, metallicities, the current 
and past star formation rates, and the fraction of light in var- 
ious structural components. Each of these features reveals 
potentially important clues for how galaxies, including their 
gas content and dark matter halos, were created and have 
evolved. These properties are also interrelated through well- 
known scalings, such as the colour-magnitude relation (e.g., 
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Sandage & Visvanathan 1978), the size-magnitude relation 
(e.g., Romanishin 1986), the Tully-Fisher relation (Tully & 
Fisher 1977), and the Fundamental Plane (Djorgovski & 
Davis 1987; Dressier et al. 1987). These correlations all de- 
scribe the scaling of structural features of normal galaxies 
with each other, or with the internal kinematics of galaxies. 

It is however not yet clear what the most basic evolu- 
tionary processes and galaxy observables are, and whether 
or not well known relations between galaxy parameters cap- 
ture essential galaxy features. Is there a single, or a group of 
processes responsible for the diversity of galaxies? Although 
it is widely known that physical properties vary along the 
Hubble sequence, galaxies within a Hubble type have a large 
range of sizes, colours, surface brightnesses, masses, and so 
on. Perhaps there is an alternative way to classify galax- 
ies based on the properties responsible for driving evolu- 
tion. One goal of this paper is determining what physical 
features of galaxies are the best for classification purposes. 
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Currently it is assumed by many that the optical morphol- 
ogy of a galaxy provides this classification, either in the form 
of a Hubble type, or through quantitative methods. However 
there is no a priori reason that optical morphologies should 
reveal information needed to classify galaxies into a mean- 
ingful system. 

We utilise fourteen properties of 22,121 galaxies in the 
nearby universe to address this problem. We begin our anal- 
ysis with the hypothesis that our fourteen measurable fea- 
tures provide the basic physical information needed to fully 
characterise nearby galaxies. We analyse these measurables 
through statistical methods, such as Spearman rank and 
principal component analyses, to determine which observ- 
ables are the most fundamental for classifying galaxies. We 
demonstrate that at a fundamental level many measured 
properties of galaxies correlate with each other. For example, 
scale features correlate strongly, and Hubble types correlate 
with many properties, including colour and stellar mass. We 
argue that this correlation is the reason that Hubble types 
have been successful in describing the nearby galaxy popu- 
lation. 

We further show through principal component analyses 
that most of the variation in galaxy properties can be ac- 
counted for by three eigenvectors dominated by scale (e.g., 
mass), the presence of star formation, and the degree of re- 
cent interaction/merger activity. As a result, our main con- 
clusion is that these three features are needed to under- 
stand the current and past evolutionary state of a galaxy. 
We further argue that galaxy mass, star formation, and 
mergers/interactions are responsible for the diversity in the 
galaxy population in the nearby universe, and should be 
used to construct a galaxy classification system. 

This paper is a more general companion, and in some 
ways a prequel, to the quantitative CAS (concentration, 
asymmetry, dumpiness) structural analysis paper presented 
in Conselice (2003). In Conselice (2003) we created a new 
quantitative way to classify galaxies under the assumption 
that mass, star formation, and recent mergers were the criti- 
cal properties for understanding galaxies. We now show that 
this assumption is justifiable through our statistical analysis 
of > 22, 000 nearby galaxies. Conselice (2003) furthermore 
developed a quantitative and reproducible method for re- 
trieving this information through the optical light of galax- 
ies. 

The outline of this paper is as follows. In §2 we give 
a summary of the physical features, and formation histo- 
ries of various types of galaxies we study. §3 describes our 
galaxy sample, and the properties we use to characterise 
each system. §4 gives a detailed description of the nearby 
galaxy population, including the relative numbers and num- 
ber densities of various morphological types. §5 presents our 
statistical analysis of our sample, while §6 utilises this infor- 
mation to produce a new physical classification system, and 
§7 is a summary. For distance measurements we use a Hub- 
ble constant of H = 70 km s _1 Mpc -1 , unless otherwise 
noted. 



2 NEARBY GALAXY TYPES AND THEIR 
POSSIBLE FORMATION HISTORIES 

There are a few generally agreed upon major nearby ba- 
sic galaxy types, generally classified according to morphol- 
ogy, whose properties and range of possible formation his- 
tories are becoming well understood. These types are: ellip- 
ticals, spirals, lenticulars, dwarf spheroidals/ellipticals, low 
surface-brightness galaxies, irregulars, and mergers, which 
we review below. These groups can, but do not necessarily 
include, other types of galaxies such as starbursting systems 
(e.g., ultra luminous infrared galaxies, ULIRGs), galaxies 
with active galactic nuclei (quasars, radio galaxies, etc.). 

Most ellipticals, particularly those in clusters, contain 
old stellar populations that likely formed early in the uni- 
verse (Kuntschner & Davies 1998). Many of these galax- 
ies are thought to form by the mergers of smaller galaxies 
(Barnes & Hernquist 1992), or from a rapid collapse of non- 
rotating gas (van Albada 1982). The merger idea for the ori- 
gin of ellipticals is currently popular, with both indirect and 
direct supporting evidence (e.g., Schweizer & Seiter 1988; 
Conselice 2006). 

Spiral galaxies are a mixed galaxy class, but all gen- 
erally show evidence of recent star formation and past 
merger/accretion activity (e.g., Kennicutt 1998; Zaritsky 
& Rix 1997). The stellar populations in spirals are also 
younger, on average, than those found in giant ellipticals 
(e.g, Trager et al. 2000). 

The lenticular class (SOs) are disk-like galaxies with- 
out the presence of strong spirals arms, star- formation, ob- 
vious dust, or gas usually associated with spirals and other 
star-forming galaxies. SO galaxies are usually found in galaxy 
clusters (Dressier 1984), but can also exist in the field (e.g,. 
van den Bergh 1997). It has been suggested since Spitzer & 
Baade (1951) that SOs originate from spirals due to environ- 
mental effects. Evidence for this includes spirals in nearby 
clusters that are highly depleted in gas (Haynes & Giovanelli 
1986), and might be evolving into SO galaxies. 

Dwarf ellipticals (dE), the most common type of 
galaxy in the nearby universe, are almost always located in 
dense environments such as galaxy clusters and groups (Fer- 
guson & Binggeli 1994). These objects resemble ellipticals 
due to their symmetric structure, and elliptical isophotes, 
although they differ in several significant properties, such 
as containing different light profiles, and how they follow 
scaling relations (e.g., Wirth & Gallagher 1984; Graham 
& Guzman 2003). These objects likely have several forma- 
tion mechanisms, including a primordial origin and evolution 
from stripped galaxies (Conselice et al. 2001,2003c). 

Low surface brightness galaxies (LSBs) are one 

of the last major class of galaxies discovered to date (Dis- 
ney 1976). These objects come in many shapes and sizes, 
but generally mimic the Hubble sequence in terms of bulge 
to disk ratios, but at a much lower surface brightness, at 
/i > 24 mag arcsec -2 (see Impey & Bothun 1997). It has 
been suggested from several pieces of observational evidence 
that LSB galaxies have had few, if any, major merging or 
interaction events with other galaxies, and are found in rel- 
ative isolation (e.g., Knezek 1993; Mo, McGaugh & Bothun 
1994). LSB galaxies therefore potentially represent a class 
of objects formed mostly through the collapse of gas, with 



The Fundamental Properties of Galaxies and a New Galaxy Classification System 3 



gradual and continual ongoing star formation, with few in- 
teractions throughout their history. 

Irregulars are typically smallish star forming galaxies 
that often lack a distinct spiral structure. For our purposes, 
irregulars include dwarf irregulars, blue compact dwarfs, and 
classical irregulars such as NGC 4449. Star formation domi- 
nates the appearance and evolution of these small, but com- 
mon, galaxies that have a wide range of surface brightnesses 
(see Gallagher & Hunter 1984 for a review). They are also 
the bluest type of non-interacting galaxy known (Huchra 
1977). A large amount of attention has been spent studying 
the dlrr members of the local group, although these objects 
also exist in clusters (e.g., Binggeli, Sandage & Tammann 
1985; Gallagher & Hunter 1986). 

Finally, interacting and merging galaxies, while 
rare in the nearby universe, were much more common in 
the distant past when roughly half of all massive galax- 
ies were undergoing major mergers (e.g., Conselice et al. 
2003). Mergers are likely a critically important galaxy for- 
mation/evolution process especially if the first galaxies were 
all of low mass. Mergers are also observationally more com- 
mon at earlier times (e.g., Conselice et al. 2003a, Bundy et 
al. 2004, Lin et al. 2004, Pope et al. 2005). 

3 NEARBY GALAXY ANALYSIS 
3.1 RC3 Sample and Data 

The sample of galaxies with which we performed our statisti- 
cal analyses were taken from galaxies in the Third Reference 
Catalog of Bright Galaxies (hereafter RC3; de Vaucouleurs 
et al. 1991). These RC3 galaxies were chosen based on their 
size (> 1'), apparent magnitude Bt < 15.5 and radial veloc- 
ity (v < 15,000 km s _1 ), although the catalogue is supple- 
mented by objects that met only one of these criteria. This 
catalog is however inhomogeneous, which we deal with in 
several ways, described in the following analysis. For each 
galaxy we culled from the literature, or derive, the following 
features. Our total sample consists of 22,121 galaxies dis- 
tributed across the sky as shown in Figure 1. The average 
galaxy in our sample is at L„, and its average galaxy colour 
is (B - V) = 0.66. 

I. T-Type (T): The mean de Vaucouleurs (1959) T-type 
is an expression of the gross morphology of a galaxy. The 
T-type values we used are means for each galaxy, com- 
puted from a variety of sources (RC3). These T- types can be 
equated with Hubble types as: -6 through -4 are ellipticals, 
-3 through -1 are SOs, and (0,1,2,3,4,5,6,7,8,9) are (SO/a, Sa, 
Sab, Sb, Sbc, Sc, Scd, Sd, Sdm, Sm), respectively. The RC3 
is the largest eye-ball estimated morphological catalogue to 
date. 

II. Size (R): The sizes of each galaxy is computed from 
the isophotal major axis diameter (D25) measured in arc- 
minutes and converted into kpc by using the distance mod- 
ulus for each galaxy. While Petrosian radii, or a radius which 
does not use isophotes, is more ideal, the effects of using an 
isophotal radii for these nearby systems is minor. 

III. Mean Effective Surface Brightness (< ji > e ): The 
mean effective surface brightness is the surface brightness 
within the half light radius of each galaxy. This parameter, 
effectively combines the size and magnitude of a galaxy and 
gives us a rough idea of the projected density of stars. 
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Figure 1. The spatial distribution of the galaxies used for our 
statistical study of galaxy populations and their properties. The 
U-shaped gap is the galactic plane, and the structure at the centre 
constitutes the Virgo superclustcr. 



IV. (U — B)o and (B — V)o colours: these are photomet- 
rically determined colour indices corrected for galactic and 
internal extinction. We used these indices as a measure of 
the average luminosity weighted age of stars in galaxies, al- 
though for older (> 6 Gyr) stellar populations these colours 
can be dominated by metallicity effects (e.g., Worthey 1994). 

V. HI linewidths (W 20 and W50): The HI line widths 
at the 20% and 50% of the maximum level in the 21-cm line 
profile. 

VI. HI linewidth ratio (J = W20/W50). We used this 
number, based on the arguments in Conselice, Bershady & 
Gallagher (2000b) as a measure of the importance of recent 
interactions and mergers. Values of I> 1.5 are interpreted as 
systems that have recently undergone interactions or merger 
with other galaxies. We discuss this index in more detail in 
§3.2 and Conselice et al. (2000b). 

VII. Central velocity dispersion (Vs). This is usually 
only defined for early type galaxies, and Vg times a size can 
be used as a measure of early type galaxy mass. 

VIII. HI line index (HI) is defined as rri2i = 16.6 - 2.5 
log Sh , where Sh is the flux of the HI line measured in units 
of W m -2 . We further subtract this value by the distance 
modulus (DM) to obtain physical quantities that are inde- 
pendent of distance. We use the HI line index as a measure 
of the cold neutral gas content in a galaxy. 

IX. Far infrared index (FIR) is defined in a similar man- 
ner as the HI line index (HI), as mpiR = —20 — 2.5 log(FIR) 
-DM, where FIR = 1.26[2.58 x/ I/ (60/xm) + / v (100jum)] 1 and 
fi/(60/im) and f„(100/im) are the total flux densities in the 
IRAS 60 and 100 urn bands (RC3) and where DM is the 
distance modulus of the galaxy. 

X. Luminosity class (A): A is the van den Bergh (1966) 
DDO, morphologically determined luminosity class index, 
which have been extended to fainter magnitudes in the RC3. 
The luminosity class values we used are weighted averages, 




4 Christopher J. Conselice 



just as the T-types, and have values that range from A = 1 
to 11 with higher A values for fainter galaxies. 

XI. The maximum rotational velocity (V max ): V ma x 
is defined as the maximum rotational velocity of a galaxy, 
as measured from HI or Ha rotation curve profiles. The size 
of a disk galaxy times V^ax i s sometimes used as a measure 
of its total mass. 

XII. The absolute B-band magnitude (Mb). 

XIII. The environmental galaxy density p. This is de- 
fined as the number of galaxies within 3 Mpc and < 1000 km 
s _1 of each galaxy in our catalog. To avoid severe incom- 
pleteness for faint galaxies, we only compute this number for 
bright galaxies with Mb < —20 and at declinations 8 > —30. 

XIV. Galaxy Stellar Mass. Stellar masses for our sample 
were computed by deriving the stellar (M/L) ratio in the B- 
band using the (B— V)— M/L ratio relation discussed in Bell 
et al. (2003), using their 'diet' Salpeter initial mass func- 
tion. There is a general correlation between stellar (M/L) 
and colour, such that bluer colors have lower stellar (M/L) 
ratios. We used this relation, and the measured (B — V)o 
colours for our galaxies to derive the (M/L) in the B-band 
for our sample. We then compute stellar masses by multi- 
plying the luminosity in the B-band for each galaxy by its 
derived stellar (M/L) ratio. 

Not all galaxies in our sample have each of these pa- 
rameters measured, therefore when comparing different pa- 
rameters we often had to use a subset of the data. We also 
associate with these parameters an error that reflects the 
variation in the measured values. 

3.2 The Galaxy Interaction/Merger Index (I) 

The galaxy interaction and merger index, /, defined by the 
ratio of the width of the HI line at 20% and 50% of max- 
imum (property VI in §3.1) is used in this paper as the 
principle indicator for finding galaxies undergoing mergers 
and interactions. It is a dynamical indicator and has been 
used at least once before for this purpose in Conselice et al. 
(2000b). However, it is not a standard galaxy measure and 
it is worth describing it, and the galaxies it selects, in more 
detail. 

The idea behind its use is that a normal galaxy with HI, 
will present a regular HI profile that is relatively narrow, or 
at least does not significantly change shape. If a galaxy is 
interacting or merging, and HI gas is disturbed by these 
dynamical events then the HI profile will have a shallow 
rise, or will have wings. This leads to higher values of I 
=W2o/W5o- Conselice et al. (2000b) show that the / index 
correlates well with the asymmetry parameter, such that 
galaxies with very large asymmetries due to merging have 
large / values. 

However, beyond the sample in Conselice et al. (2000b), 
this index has not been tested to determine its success rate 
in finding merging galaxies. To understand this better we ex- 
amined the literature and Digitised Sky Survey (DSS) imag- 
ing for several hundred systems with I > 2. What we found 
was a diversity of galaxy morphology, from spirals and ob- 
vious mergers, to galaxies classified as SOs/ellipticals. How- 
ever, it was nearly always the case that a system with a reg- 
ular apparent morphology was in a galaxy pair, such that 
it is likely involved in an interaction with another galaxy. 
Examples of the optical morphology of the high-/ systems 
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Figure 2. Examples of galaxies within our sample, mostly taken 
from the DSS, which have large 20% to 50% HI line width ratios, 
/. As can be seen galaxies in interactions, in pre-mergers and 
post-mergers will all produce high HI line width ratios. 

taken from the NASA/IPAC Extragalactic Database (NED) 
and the DSS are shown in Figure 2. It is clear from these 
images that the I indicator is a more sensitive probe than 
the asymmetry index (Conselice 2003), which is designed 
to find galaxies that have already merged. It appears that 
the I index will become large during an interaction as well 
as throughout the merger of two galaxies. There is however 
some contamination from non-interacting galaxies into our 
high / sample, as might be expected. This rate is however 
low, and does not affect our results, as it is only at the 5% 
level. Much of this contamination seems to arise from edge- 
on disk galaxies. 



4 NEARBY GALAXY PROPERTIES 

Our primary goal using our large sized-limited sample 
of nearby galaxies with morphologies is determining the 
strongest correlations between measurable parameters (§5). 
Before doing this, we examine the nearby galaxy population, 
its overall distribution of properties, and how each relates 
to morphological type. One purpose for this is having a well 
defined morphological census of nearby galaxies to compare 
with more distant samples studied with the Hubble Space 
Telescope. 

4.1 Galaxy Stellar Mass and Spectro-Morphology 

at z ~ 

As the Hubble Space Telescope can now reveal the morpho- 
logical properties, and stellar masses of thousands of galaxies 
at z > 0.5 (Papovich et al. 2005; Bundy et al. 2005; Con- 
selice et al. 2004, 2005a), direct morphological and stellar 
mass comparisons with the z — galaxy population are 
now possible (e.g., van den Bergh, Cohen & Crabbe 2001; 
Conselice et al. 2005). It is thus desirable to determine the 
morphological distribution of galaxies at z ~ and how 
these relate to stellar masses and colours. 

We first calculate some basic spectro-morphological 
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Table 1. Morphological Fractions and Densities of Nearby Bright Galaxies 
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Figure 3. The radial velocity distribution of the galaxies used 
in our analysis verses their absolute magnitude (Mb). The box 
formed by the solid line outlines the nearly complete sample we 
used for analyses that require near complete samples brighter 
than M B < -20. 

properties of nearby galaxies, including how stellar masses 
relate to morphology at z ~ 0. Our catalog of 22,121 RC3 
galaxies is however inhomogeneous, particularly for low sur- 
face brightness systems. As Figure 3 shows, we sample 
fainter galaxies at lower radial velocities, implying we are 
incomplete at fainter magnitudes. To partially overcome the 
incompleteness of our catalog, we limit our morphological 
statistical analyses to systems which have Mb < — 20 and 
v < 10,000 km s' 1 (z ~ 0.033, 133 Mpc) (shown by the 
lines in Figure 3). When investigating fainter galaxies with 
Mb < —18, we limits our analyses to objects with velocities 
v < 3000 km s _1 (« ~ 0.01, 40 Mpc). 

In addition to the Mb < -20 and v < 10, 000 km s _1 
cut, we furthermore only consider galaxies between declina- 
tion -50 < 8 < 50, and take into account the incompleteness 
of the galaxy distribution due to the Milky Way (Figure 1) 
when calculating total number densities. Note that these 
number densities are based on the RC3 selection, and are not 
complete in any sense. Our densities however generally agree 
well with the densities derived from the morphological-type 
specific Schechter functions published in Marzke et al. (1998) 
and Nakamura et al. (2003). We investigate the complete- 
ness of our sample down to Mb = —20 through statistical 



Figure 5. The fraction of different morphological types for all 
nearby bright galaxies with Mb < —22 (solid line), Mb < —20 
(dashed line) and M B < -18 with V < 3000 km s _1 (dotted line). 
The most common galaxy in our sample is therefore a mid-type 
spiral, which docs not change significantly at different luminosity 
cuts. 

tests. From this sample we create V/Vmax diagrams for all 
morphological types, revealing that we are nearly complete 
to 10,000 km s _1 for galaxies brighter than Mb = —20. 

We show fractional histograms in Figure 4 of the dis- 
tribution of morphological types, (B — V)o colours, stellar 
masses, and absolute magnitudes (Mb) for the entire sample. 
Figure 5 shows the fraction of all galaxies that are classified 
into different T-types. From Figure 5, the most common 
galaxy type in the nearby universe appears to be mid/late 
type spirals, with a peak fraction at T-type = 4 (Sbc) . This 
peak remains unchanged when we cut the sample at brighter, 
M s < -22, or fainter, M s < -18, limits. 

The relative percentages of nearby galaxies with differ- 
ent morphological types are listed in Table 1. It appears 
that ~ 75 % of all nearby bright galaxies with Mb < —20 
are spirals, with ~ 40% late type spirals (Sb-Sc). Thus, in 
the nearby universe most bright galaxies (with Mb < —20) 
are spirals with small bulges. After accounting for spatial 
incompleteness in the catalog, we further find that the total 
number density of spirals is ~ 3 x 10 6 Gpc~ 3 at Mb < —20 
(Table 1). 

Early-type galaxies, including ellipticals and SOs, ac- 
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Figure 4. The distribution of morphology (T-type), colour, absolute magnitude, and stellar mass (M,) for all galaxies used in this study. 
We rapidly become incomplete at Mb > —20 and at M < 5 X 10 10 Mq within our full volume (see Figure 3). As can be seen, we span 
the entire range of Hubble types and galaxy colours. 



count for nearly all of the remainder (22%) of nearby bright 
galaxies, with the rest ~ 2% irregulars (Table 1). These num- 
ber densities and morphological type fractions are being di- 
rectly compared to bright galaxies at high redshifts where 
these same properties can now be determined (e.g., Con- 
selice 2004; Conselice et al. 2005a) . Figure 6 shows the frac- 
tion of each morphological classed type that is also classified 
as having a bar, ring or part of a multiple system, which we 
discuss below. 

Another feature we investigate is how galaxy morphol- 
ogy correlates with the positions of galaxies in the colour- 
magnitude plane (Figure 7). There is a well characterised 
red-sequence and blue-cloud in the (U — B)o vs. Mb dia- 
gram (e.g., Baldry et al. 2004), which has been proposed as 
a fundamental feature for understanding how galaxies have 
evolved. When we label the morphological types of galaxies 
on the colour- magnitude digram we find a good correlation 
between the morphologies of galaxies and their positions. 
The early types (denoted by solid red squares) are mostly 
on the red-sequence, while later types (T-type > 3, blue 
crosses) are in the blue cloud, and mid-types (T-type > & 
T-type < 3) are in between. 

This is also the case for stellar masses, with the M* > 
10 11 M Q galaxies mostly ellipticals on the red-sequence. 
Lower mass galaxies are increasingly bluer. This can also be 
seen in Figure 8, where we plot the morphological fraction as 
a function of stellar mass. At the highest masses, nearly all 
galaxies are ellipticals, but this gradually changes, and by 
M* = 10 11 M ellipticals and spirals each represent about 



~50% of the population. At the low mass end, irregulars be- 
gin to dominate the population, with spirals accounting for 
only 20% of the population at M* < 10 9 M Q and ellipticals 
accounting for about 10%. 



4.1.1 Bars 

Bars and rings are usually, but not uniquely, associated with 
spiral galaxies (e.g., Buta & Combes 1996). Bars are dynam- 
ical features that reveal, and contribute to, several galaxy 
evolution processes such as star formation and AGN activ- 
ity (Barnes & Hernquist 1992), or galaxy interactions. It is 
debated however whether bars are a fundamental feature of 
galaxies that should be used in classifications (e.g., Abra- 
ham & Merrifield 2000; Conselice 2003). The evolution of 
bars through cosmic time is also debated, although there 
are some claims that the bar fraction and sizes of disks have 
not evolved since z ~ 1 (Jogee et al. 2004; Ravindranath et 
al. 2004). As the evolution of bars will continue to be ad- 
dressed in the future through observations with the Hubble 
Space Telescope, we compute basic physical information on 
barred galaxies in the nearby universe for comparison pur- 
poses (Figure 6a, Table 2). 

Table 2 lists, for all barred galaxies in the nearby uni- 
verse, the fractional breakdown of host galaxy morpholog- 
ical type, and the number densities of each type, for sys- 
tems brighter than Mb < —20. For example, 43% of nearby 
bright barred galaxies are found in late type disks, and late 
type disks with bars have a total number density of 5.2±1.4 
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Table 2. Bar, Ring and Multiple Component Fractions as a Function of Morphology for Galaxies with Mb < —20 



Bar Ring Multiple 

Type (# Gpc -3 hf Q0 ) Fraction Type (# Gpc" 3 h 3 00 ) Fraction Type (# Gpc" 3 h 3 00 ) Fraction 



E 


1.2±0.5 xlO 4 


0.01 


E 


9.6±10 xlO 2 


0.00 


E 


3.0±1.1 XlO 4 


0.11 


SO 
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5.2±1.4 xlO 5 


0.43 


IDisk 


1.4±0.4 xlO 5 


0.37 


IDisk 


9.3±3.0 xlO 4 


0.36 


Irr 


2.3±0.9 xlO 4 


0.02 
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4.7±3.0 xlO 3 


0.01 
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1.2±0.6 xlO 4 
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Figure 6. The fraction of different galaxy types that have morphological evidence for the presence of bars, rings, and multiple components. 
These fractions are also plotted as a function of luminosity, with Mb < —22 (solid line), Mb < —20 (long-dashed line) and Mb < —18 
with V < 3000 km s" 1 (dotted line). As can be seen, the bar fraction is roughly constant at T-types > 0, while ringed galaxies are 
typically early type spirals, and galaxies with a multiple component classification tend to span all Hubble types. 



xlO 5 hfoo galaxies Gpc -3 . Not surprisingly, we find that al- 
most all (97%) of barred galaxies are classified as disks/SOs, 
with an equal proportion belonging to early and late types. 
Roughly 10% of all barred galaxies are SOs, while 3% of 
classified barred galaxies are ellipticals and irregular galax- 
ies. The small fraction of bars that are found in ellipticals 
shows that these systems have been misclassified, as a bar 
needs a disk to exist. 

By examining the fraction of galaxies in each morpho- 
logical type that are barred (plotted in Figure 6), we find 
that roughly 50% of all disk galaxies are barred, while the 
fraction of barred galaxies increases for later types. There is 
also a fairly high fraction of nearby galaxies with bars, with 
a number density in the nearby universe of ~1.2 xlO 6 hf 00 
barred galaxies Gpc -3 . 

4.1.2 Rings 

Rings often follow the presence of bars in galaxies (Buta & 
Combes 1996), and this can be quantitatively demonstrated 
as they occupy the same types of galaxies. Of all ringed 
galaxies brighter than Mb < —20, nearly all (98%) are lo- 
cated in disks/SOs (Table 2). The fraction of ringed galaxies 
which are hosted inside ellipticals, SOs and irregulars is simi- 
lar to the fraction with bars. There are however fewer ringed 
galaxies than barred galaxies in the nearby universe, with 
a number density of 3.8 xlO 5 ringed galaxies Gpc -3 . Fig- 
ure 6b shows the fraction of galaxies of all morphological 



types with rings, revealing that early disks are more likely 
to host a ring than later types. 



4-1.3 Multiple Components 

In addition to bars and rings, the RC3 also includes a clas- 
sification for galaxies that have multiple components. Intu- 
itively, this might correlate with the presence of merging 
activity. About 71% of galaxies classified as having multiple 
components are disks, while 11%, 13% and 5% of multiple 
galaxies are ellipticals, SOs and irregulars, respectively. How- 
ever, a different pattern emerges when we consider the frac- 
tions of each galaxy type that have a multiple classification. 
Figure 6c plots the fraction of galaxies in each T-type that 
have morphological evidence for multiple components. This 
demonstrates that only about 5 - 10% of all galaxies in the 
RC3 are classified as multiple galaxies. A higher fraction of 
irregulars are classified as having multiple components, al- 
though this is likely partially due to low number statistics. 
If these fractions are tracing the merging activity of galaxies 
in the nearby universe, the derived fraction of 5% is close to 
the z ~ merger fraction estimates from pairs (e.g., Patton 
et al. 2000). 
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Table 3. Spearman Correlation Matrix For All Galaxies 



T-Type R < fj, > e (U-B) (B-V) W 20 V s HI FIR A V max M B I p M, 



T-type 


1.00 


-0.24 


0.53 


-0.79 


-0.80 


-0.51 


-0.43 


0.18 


0.13 


0.91 


-0.50 


0.25 


0.23 


-0.14 


-0.68 


R 


-0.24 


1.00 


0.08 


0.23 


0.20 


0.66 


0.48 


-0.76 


-0.55 


-0.62 


0.65 


-0.85 


-0.36 


-0.13 


0.51 
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0.53 


0.08 


1.00 


-0.35 


-0.36 


-0.38 


0.09 


-0.01 


0.08 


0.56 


-0.35 


0.00 


0.07 


-0.10 


-0.07 


({/ - B) 


-0.79 


0.23 


-0.35 
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0.91 


0.51 


0.60 


0.04 


0.03 


-0.67 


0.48 


-0.23 


-0.30 


0.19 


0.75 


(B - V) 


-0.80 


0.20 


-0.36 


0.91 


1.00 


0.47 


0.57 


0.09 


0.03 


-0.64 


0.47 


-0.19 


-0.28 


0.22 


0.80 


W20 


-0.51 


0.66 


-0.38 


0.51 


0.47 


1.00 


0.30 


-0.54 


-0.39 


-0.72 


0.87 


-0.68 


-0.38 


-0.03 


0.51 


V S 


-0.43 


0.48 


0.09 


0.60 


0.57 


0.30 


1.00 


0.09 


0.07 


-0.31 


0.21 


-0.55 


-0.16 


0.06 


0.67 


HI 


0.18 


-0.76 


-0.01 


0.04 


0.09 


-0.54 


0.09 


1.00 


0.56 


0.53 


-0.55 


0.68 


0.23 


0.2 


-0.05 


FIR 


0.13 


-0.55 


0.08 


0.03 


0.03 


-0.39 


0.07 


0.56 


1.00 


0.40 


-0.46 


0.64 


0.05 


0.23 


-0.08 


A 


0.91 


-0.62 


0.56 


-0.67 


-0.64 


-0.72 


-0.31 


0.53 


0.40 


1.00 


-0.72 


0.67 


0.40 


-0.03 


-0.51 


Vmax 


-0.50 


0.65 


-0.35 


0.48 


0.47 


0.87 


0.21 


-0.55 


-0.46 


-0.72 


1.00 


-0.69 


-0.36 


-0.05 


0.55 


M B 


0.25 


-0.85 


0.00 


-0.23 


-0.19 


-0.68 


-0.55 


0.68 


0.64 


0.67 


-0.69 


1.00 


0.31 


0.10 


-0.59 


I 


0.23 


-0.36 


0.07 


-0.30 


-0.28 


-0.38 


-0.16 


0.23 


0.05 


0.40 


-0.36 


0.31 


1.00 


-0.01 


-0.28 


P 


-0.14 


-0.13 


-0.10 


0.19 


0.22 


-0.03 


0.06 


0.20 


0.23 


-0.03 


-0.05 


0.10 


-0.01 


1.00 


0.13 


M» 


-0.68 


0.51 


-0.07 


0.75 


0.80 


0.51 


0.67 


-0.05 


-0.08 


-0.51 


0.55 


-0.59 


-0.28 


0.13 


1.00 



Spearman correlation matrix for all the galaxies in our sample. Displayed are the correlation between various physical properties, including 
morphology (T-type), radius (R), surface brightness (< fi > e ), colour (both (U — B)q and (B — V)o), HI line width (W20): central velocity 
dispersion (Vs), HI and far infrared line fluxes (FIR), the luminosity class (A), the maximum rotation velocity (V m ax), the absolute magnitude 
(Mb), an interaction index (/), local density (p) and the stellar mass M*. The numbers in bold show the strongest correlations between the 
various parameters. Tables 4-6 show these correlations after the same is divided into early/late, red/blue and bright/faint systems. 




Figure 7. The colour-magnitude relation for galaxies in our sam- 
ple. The solid line divides red-sequence galaxies from blue cloud 
galaxies (Faber et al. 2005). Early types are denoted by solid 
red squares, while later types (T-type > 3) are labelled as blue 
crosses, and mid- types (T-type > &; T-type < 3) are labelled 
as green triangles. 



5 GALAXY PROPERTIES AND THEIR 
CORRELATIONS 

5.1 Statistical Analysis of the Nearby Galaxy 
Population 

We now examine how the 14 properties described in §3.1 
correlate with each other, and how strongly. We do this 
by performing Spearman rank correlation tests on various 



Figure 8. The fractional contribution of various morphological 
types to the local stellar mass density down to a lower limit of 
M„ = 10 s Mq . The solid line shows the ellipticals, the dashed 
line spirals, and the dotted line irregular galaxies. 



cuts of our data (Tables 3-6). These correlation analyses 
are performed on the entire sample (Table 3), when di- 
vided by early and late morphological types (Table 4) , colour 
at (B — V)o = 0.7 (Table 5), and absolute magnitude at 
Ms = —20 (Table 6). Our sample for the Spearman rank 
correlations include the RC3 sample out to radial velocities 
of 15,000 km s _1 . We discuss in specific cases where not hav- 
ing a volume limited sample may be affecting our results. 

Our entire sample's correlation matrix (Table 3) reveals 
several strong relationships, including between: Vs and Mb 
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(Faber- Jackson), log V max and Mb (Tully-Fisher) , and size 
(R) with Mb (Kormendy relation). However, the relations 
we find are weaker than these scaling relationships, as the 
sample here has not been divided into different morpholog- 
ical types. In general however, there are many correlations 
between the 14 parameters, revealing that to one degree 
or another most properties of galaxies are interrelated. Is 
this interrelation the result of a single principle, or a group 
of physical principles producing different observational ef- 
fects? For example, it is possible that size and Mb have a 
strong correlation because they both fundamentally scale 
with an underlying feature, which is likely the total mass? 
The same is also true for the strong correlation between both 
Mb and size with internal velocities - as things get bigger 
and brighter they are likely to have more mass, an effect 
traced by larger internal motions. 

We divide these parameters into three classes, following 
the ideas described in Bershady, Jangren & Conselice 2000 
(hereafter BJC00). BJC00 determined from an analysis of 
nearby galaxy images that galaxies can be classified accord- 
ing to scale, form, and spectral index. In this study colour 
represents the spectral index, scale includes the size, vari- 
ous fluxes, internal velocities and stellar mass, while form is 
given by the T-type, bars and the multiple component clas- 
sification. We investigate which of these measurable prop- 
erties are the most fundamental, and which are the most 
useful for classifying galaxies. 

5.1.1 Form Correlations 

We use the T-type to represent the form of a galaxy and to 
first order assume it correlates with the bulge to disk ratio 
of each galaxy (see also Conselice 2003). T- types (or Hub- 
ble types) have been known since at least Holmberg (1958) 
to correlate broadly with other physical properties, such as 
recent star formation. Almost by definition this should be 
the case, as ellipticals and SOs, which populate the earlier T- 
types, appear fundamentally different from later type galax- 
ies. This difference has been described in more detail by e.g., 
Roberts & Haynes (1994). It is not yet clear however which 
physical properties of galaxies correlate strongest with mor- 
phology, and if Hubble morphology is enough to describe 
the most fundamental aspects of galaxies. With our Spear- 
man correlation coefficients, we can begin to address this 
question. 

Table 3 shows the Spearman coefficients for the entire 
RC3 galaxy sample. If we consider strong correlations those 
with coefficients > |0.5| (listed in bold) then T-types corre- 
late strongly with both (B — V)o and (U — B)o colours, W20, 
the luminosity index (A), V max , Mm, and M„. This reveals 
several things about what a T-type is measuring. Since T- 
types correlate with B/T ratios, then the strong correlation 
between T-types and colour is not surprising, as it is well 
known that the relative contributions of bulge and disk light 
correlates with old and young stellar populations. The cor- 
relation between the stellar masses of galaxies and Hubble 
T-type (-0.68 coefficient), is stronger than the correlation be- 
tween T-type and absolute magnitude (0.25), showing that 
fundamentally the Hubble sequence (or bulge to disk ratio) 
is partially one of stellar mass. There is likewise a strong 
correlation (coefficient 0.75-0.8) between stellar mass and 
colour (Table 1), which we discuss in §5.1.2. 



T-types also correlate strongly with other scale fea- 
tures, as the strong correlations between T-types and V max , 
W20 , Vs and M* demonstrates. The absolute magnitude and 
size of galaxies do not show the same strong correlation, 
although later type galaxies are smaller and fainter. This 
weaker correlation between T-type and magnitude/radius 
might be the result of magnitude/radii being secondary scale 
features. This is because the measured size of a galaxy is 
influenced by both the form of a galaxy (i.e., how light is 
distributed, for example bulge vs. disk components) and the 
intrinsic scale, and magnitude is influenced by scale, and 
spectral type, as younger, bluer stars will make a galaxy 
brighter. In this sense these values are not clean measure- 
ments of the scale. However, the internal velocities of galax- 
ies and stellar masses are not affected by these secondary 
effects, and they thus represent more accurately the scale of 
a galaxy. 

Table 4 shows the Spearman correlation matrix between 
our 14 parameters when the sample has been divided into 
early (T < 1) and late type (T > 1) galaxies. We repeat 
the Spearman correlation analysis on both early and late 
types, to determine how correlations between T-types and 
the other parameters are affected by this division. Table 4 
shows that for early-types galaxies with T < 1 there is very 
little to no correlation between T-type and other physical 
parameters, with the exception of perhaps colour. This is 
such that later-type early-types (E6-7,S0) are bluer and have 
a slightly lower stellar mass. 

When we examine correlations between T-types and 
physical properties for the late-type galaxies (Table 4) many 
of the correlations seen in Table 3 remain, but with most 
weaker. The correlation with colour for the late types de- 
clines from ~ —0.8 to —0.61, as does the correlation with 
stellar mass. Interestingly, the mild correlation between T- 
types and surface brightness (fi E ) seems to remain for the 
late types, such that later-types have a lower surface bright- 
ness. This is perhaps the only correlation in the Hubble se- 
quence (besides Vmax) which is driven mostly by later type 
galaxies alone. 

Most correlations between physical features and Hub- 
ble types is strongest when considering the entire sequence. 
One exception to this is the correlation between stellar mass 
and colour which remains strong after dividing the sample 
into early/late Hubble types (Table 4). This shows that this 
correlation is not driven entirely by morphology. There is 
however little correlation within the early-types themselves, 
with the exception of scale-scale correlations. This demon- 
strates that the division of early types onto the Hubble se- 
quence, which is typically done by estimating the axis ratios 
of galaxies (Sandage 1961), is largely devoid of much useful 
physical information (Kormendy & Bender 1996). The fact 
that T-types correlate strongly with colour and stellar mass 
can be understood to first order if star formation follows the 
presence of disks. Galaxies with types T > 1 have a diversity 
in B/T ratios, and a mix of segregated stellar populations 
(Papovich et al. 2003), while ellipticals have B/T ~ 1 and 
tend to have nearly homogenous old stellar populations. 

5.1.2 Spectral Type Correlations 

As discussed in §5.1.1 galaxy T-types correlate strongly with 
colour. Spectral type, as denoted by either the (B — V)o 
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Table 4. Spearman Correlation Matrix For Early/(Late) Galaxies 



T-type 1. 

-0 

R -0. 

-0. 

</i>e -0 

-0 

(U - B) -0. 

0. 

(B - V) -0. 

0. 

W20 O.i 
-0. 

v s -0. 

0.' 

HI -0 
0. 

FIR -0 
1. 

A 0. 

0. 

Vmax 0.' 

-0. 

M B 0. 

0. 

/ -0 
0. 

p -0. 

0. 

M» -0. 

-0. 

or (U — B)o colour, also correlates well with the internal 
velocities of galaxies, i.e. W20, V max , and Vs as well as stellar 
masses (Table 3). The colours of galaxies however do not 
strongly correlate with any other galaxy property, including 
sizes, HI/FIR luminosities and absolute magnitude. 

Late-types are blue, and systems with lower internal 
velocities, and lower stellar masses are also bluer, as dis- 
cussed above. There is no a priori reason for galaxies with 
lower internal velocities, lower stellar masses, to have av- 
erage younger stellar populations, while more massive sys- 
tems have a larger fraction of older stars. This correlation 
between ongoing star formation, as revealed through bluer 
colours, and stellar mass is a form of the downsizing seen 
up to z ~ 1 (Bundy et al. 2005, 2006), which is still an 
unexplained physical effect that dominates the way galaxies 
formed. 

We separated the RC3 sample into red and blue sys- 
tems, divided at (B — V)o = 0.7, and then redetermine 
how these various correlations change (Table 5). We find 
a similar effect as when we divide the sample by morphol- 
ogy. Red galaxies show fewer correlations between different 
parameters (Table 5), again with the exception of colour, 
such that later-types with red colours, are bluer. However, 
as for early-types, there is little correlation between colour 
and stellar mass for red galaxies. Star forming blue galaxies 



.29) 
0.78) 
L09) 
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1.23) 
i.20) 
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have physical properties that correlate well with each other, 
i.e. stellar mass with colour, HI with stellar mass, and size 
with Vmax- Interestingly, the correlation between T- types 
and colour within the red galaxies is slightly stronger than 
it is for the blue galaxies (~ —0.55 vs. —0.45). There is 
still a fairly strong (~ 0.5) correlation between (B — V)o 
colour with internal velocity dispersion (Vs), and size with 
Vs (0.56), for red galaxies, revealing that these relationships 
are independent of star formation. 

5.1.3 Scale 

As was mentioned earlier, the different scale features of 
galaxies, namely: size, stellar mass, internal velocities, mag- 
nitude, HI and FIR magnitudes, all correlate with each other 
with coefficients ~ 0.5 — 0.7 for most correlations. Some of 
these strong correlations are: R with: W20 (0.66), Vs (0.48), 
V max (0.65), M B (-0.85), and M, (0.51), and M B with: W20 
(-0.68), Vs (-0.55) and V max (-0.69). This indicates that 
measures of scale correlate with each other, and are inde- 
pendent of classification. There is also a correlation between 
scale parameters, particularly internal velocities, and the 
morphological T-type, as discussed in §5.1.1. 

Correlations between scale parameters remain when ex- 
amining galaxies divided by morphology, colour and abso- 
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Table 5. Spearman Correlation Matrix For Red/(Blue) Galaxies Divided at (B — V)o = 0.7 



T-type 1. 

-0 

R 0.' 

-0. 

<fi> e 0. 

-0 

(U - B) -0. 
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(B - V) -0. 

0.. 

W20 O.i 
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1. 
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-0 

Vmax 0. 
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M B 0. 

0. 
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0. 

p 0. 

0. 

M» -0. 

-0. 

lute magnitude (Tables 4 — 6). For example, the correlation 
between M B and R remains quite strong no matter how the 
sample is divided with coefficients of ~ 0.85. Other strong 
correlations are between internal velocities, Vs, for early- 
types, V max for late-types and absolute magnitude (Mb), 
with values —0.68, —0.72, respectively. 

Correlations for systems divided into blue and red 
galaxies at (B— V)= 0.7, and for systems divided into 
'bright' galaxies with Ms > —20, and faint galaxies with 
Ms < —20 show a similar dichotomy as late/early types 
(Table 5-6). Likewise, as for red and early-types, the bright 
systems tend to have weaker correlations compared to when 
the entire magnitude range is considered. 

For the more luminous galaxies, parameters that corre- 
late strongly are: magnitude/stellar mass with size scalings 
and internal velocity (Vs), as well as M», colour and T-type 
among each other (Table 6). This is also the case for fainter 
galaxies where there is a strong correlation (~ 0.8) between 
M* and colour, such that more massive galaxies are redder. 
Both W20 and V m ax correlate well with most parameters 
except with the / index, an indication that for fainter sys- 
tems V m ax is measuring something fundamental about the 
scale of these systems. This is also true when examining the 
correlation of internal velocities (Vs) of the more massive 
systems. 
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0.68) 

5.1.4 Local Density 

The local density (hereafter environment) of a galaxy has 
been known since Hubble (1926) to correlate with galaxy 
populations. This morphology-density relation can be sum- 
marised by stating that in locally dense environments a large 
fraction of galaxies are early- types (either SOs or ellipticals), 
while in lower density environments spiral galaxies are more 
common (e.g., Dressier 1980; Postman & Geller 1984). 

We can reexamine the correlation of environment with 
other physical properties by using the calculated galaxy 
number density (p) parameter (§3.1). We however do not 
find a strong correlation between p and other parameters, 
including morphological type. This is perhaps because galax- 
ies of all types are found in all environments, but their 
relative proportions can change, creating for example, the 
morphology-density relation. Table 3 shows that the corre- 
lations between local density and properties such as colour 
and T-type have the general trend expected, such as later 
types and bluer galaxies are found in lower density envi- 
ronments. However, these correlations are not as strong as 
many others. 
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Table 6. Spearman Correlation Matrix For Bright / (Faint) Galaxies Divided at Mb 
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5.1.5 Summary of Correlations 

It is clear from the above discussion that there are strong 
correlations between various galaxies properties. A few of 
these correlations, such as size and brightness are very 
strong, and are independent of other properties such as re- 
cent star formation. Perhaps the most interesting correlation 
is the one between T-type, stellar mass and colour, which 
all seem to relate to each other. 

Scale-features, such as size and stellar mass correlate 
strongly, which does not change after dividing the popula- 
tion into early/late, red/blue or bright/faint systems. How- 
ever, other correlations break down when we examine them 
within these subsets, showing the scale is a fundamental fea- 
ture of galaxies. 

In particular, many correlations remain for late-type, or 
star forming galaxies, as defined by (B— V) < 0.7, but do not 
appear as strong as for the early/red galaxies. These corre- 
lations, which break down in the absence of recent star for- 
mation are: stellar mass/magnitude with colour and HI/FIR 
magnitudes, and T-type with surface brightness. There is 
also little correlation between T-types and size when galax- 
ies are broken into red/blue systems 

Ongoing star formation is therefore a critical component 
for trying to understand the properties of galaxies. Finally, 
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if we accept the idea that the strongest correlations with 
Hubble type reveal its underlying meaning, then the fact 
that T-types correlate with colour and internal velocities 
demonstrates that the Hubble sequence is fundamentally one 
of decreased total mass and increased dominance of young 
stars, which accounts for its success in the nearby universe. 

Note, that the form index, /, which has been shown 
to correlate with dynamical disturbances (Conselice et al. 
2000b), does not correlate strongly with any of the other ob- 
servational properties we study in this paper (Table 3) . This 
reveals that I is not directly related to, or always causes, 
other parameters to change, with the exception of perhaps 
the luminosity index. 

5.2 Principal Component Analysis 

We carry out a principal component analysis of the proper- 
ties of our sample to determine the most fundamental fea- 
tures needed to describe nearby galaxies. The question we 
address is whether we can use only morphological type to 
describe a galaxy's major properties. We perform several 
principal component analyses (PCA) to determine which 
features cannot be accounted for by other properties, and 
to reduce the dimensionality of our data set with its 14 pa- 
rameters. 
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Table 7. PCA Eigenvector (EV) Projections onto Galaxy Prop- 
erties 
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PCA methods are commonly used in astronomy, and 
good introductions to the topic can be found in e.g., Faber 
(1973), Whitmore (1984), Han (1995), Brotherton (1996) 
and Madgwick et al. (2003). The basic goal of a PCA analy- 
sis is to reduce the number of variables (n) to a set of m < n 
parameter combinations. This is done through the calcula- 
tion of eigenvectors and eigenvalues for the different param- 
eter combinations. The first eigenvector is a line through 
the cloud of n data points which minimises the variance. 
Remaining eigenvectors are calculated until all the variance 
is explained. Eigenvalues reveal how much of the variation is 
accounted for by each eigenvector. Typically the eigenvalue 
for an eigenvector must be ^ 1 to be significant. Before we 
carry out our various PCA analyses, we linearise all vari- 
ables, subtract the mean from each variable, and normalise 
by the standard deviation. 

We preformed a large number of PCA analyses on our 
data. A typical analysis that includes most of our proper- 
ties is shown in Table 7. We have created a new index here 
called "Vel" which is the internal velocity of our galaxies, 
Vmax, or Vs when V max is not available. When we calcu- 
late the total masses of our galaxies using ~ R x V^ax) or 
total masses derived based on scalings derived from models 
(Conselice et al. 2005b), we find the same correlations. The 
eigenvectors listed in Table 7 show that several vectors are 
needed to reach a reasonable accounting of the total vari- 
ance in galaxy properties. The total summation of the five 
eigenvectors listed in Table 7 represent 83% of the variance 
in the properties of galaxies. Due to the errors on the mea- 
surements, this is likely the limit for meaningful eigenvec- 
tors. In fact, when we plot the eigenvalues for the analysis 
in Table 7 vs. principal component rank we find a shallow 
slope past EV-5 than before it, suggesting that these higher 
components are dominated by noise. 

The first eigenvector (EV 1; Table 7) is dominated by 
Mb, M„, size (R), and HI content, which all positively cor- 
relate (Mb was converted into a positive linear flux for this 
analysis). The first vector can therefore be thought of as 
the scale of the system. The second vector is dominated by 
the T-type, the (B — V)o colour, and surface brightness. The 
strong correlation between Hubble types and colour suggests 
that this vector measures the presence of recent star forma- 
tion in a galaxy. Vector 3 and Vector 4 are both dominated 
by the merger parameter (7) and the FIR luminosity, and as 



Table 8. EVs for PCA with T-type, Colour and M» 



Property 


EV 1 


EV 2 


EV 3 


T-type 


-0.57 


-0.13 


0.60 


(B - V) 


0.63 


-0.04 


-0.10 


I 


0.06 


0.96 


0.23 


M» 


0.52 


-0.20 


0.76 


Eigenvalues 


2.0 


1.02 


0.64 




50.8% 


25.4% 


16.1% 



such it represents the merging properties of galaxies. Vector 
5 is dominated by the internal velocities of galaxies. 

We find a similar pattern after we remove various vari- 
ables from the PCA analysis. If we do not consider FIR or 
surface brightness in the analysis, we find that the eigenvec- 
tors are always dominated by I, (B — V)o, and M* in the 
first three principal components. We find in general that the 
T-type is a significant contributor to many components af- 
ter the analysis is carried out through different cuts in the 
data. 

When we do a basic PCA analysis on our sample sepa- 
rated into only T-types, a spectral index ((B — V)o colour), 
the merger index (I), and stellar mass (M*) we obtain Ta- 
ble 8. We find that the T-type, stellar mass and colour arc 
highly correlated and dominate the first component. The / 
merger index defines the second principal component, while 
the third component is dominated by the T-type, and the 
stellar mass. This shows that the T-type is successful as a 
general galaxy classification criteria as it correlates strongly 
with colour and stellar mass, and several other properties 
(cf. Table 3) . Therefore the three most fundamental proper- 
ties for describing nearby galaxies appear to be their scale, 
recent star formation, and the presence of mergers. 

From §5.1 and Table 3 we know that the interac- 
tion/merger index, /, does not correlate strongly with any 
other galaxy property, with the exception of the luminosity 
index. The meaning of this is that galaxies which are under- 
going a merger cannot be accounted for by other measurable 
features. When we analyse only galaxies undergoing merg- 
ers, defined as objects with I > 1.5, we find that three prin- 
cipal components are necessary to describe these systems. 
In fact, for systems with M > 1.5 the correlation between 
T-type and (B — V)o colours breaks down. We conclude from 
the above analyses that our hypothesis, described in Con- 
selice (2003) that mass, star-formation and the presence of 
mergers are fundamental criteria, is a reasonable assump- 
tion. 

5.3 Galaxies in PCA Space 

In this section we examine whether the principal compo- 
nent projections for our sample are correlated, and if dif- 
ferent galaxy types separate cleanly within PCA space (see 
Madgwick et al. 2003 for examples of this using galaxy spec- 
tra) . For this analysis we examine the principal components 
shown in Table 8. Figure 9 shows the correlation between the 
projections of our galaxies onto their first three eigenvectors 
listed in Table 8, which as we have discussed, generally corre- 
late with colour (pi), interactions/mergers (P2), and stellar 
mass (P3). We furthermore label the morphological type of 
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galaxies on Figure 9, which shows that this relatively sim- 
ple PCA analysis separates nicely the major galaxy classes: 
early-types, disks, irregulars and merging galaxies. The left 
panel of Figure 9 shows a correlation between the projec- 
tion of the eigenvectors dominated by recent star formation 
(colour) and scale (stellar mass). There is however a per- 
pendicular relation as well which is occupied mostly by disk 
galaxies. 

The correlation between pi and p2 shows a similar pat- 
ten, but with only a slight correlation between the two eigen- 
vector projections. However, it is clear that interacting and 
merging galaxies, with I > 2, deviate from the nominal re- 
lation. Examples of these high-/ systems are shown in Fig- 
ure 2. There is also a correlation between the second and 
third eigenvector projections for most galaxies in our sam- 
ple. Just as in the right panel of Figure 9, galaxies with a 
high / index deviate from this relation with a larger projec- 
tion onto the second eigenvector. 



5.4 The Strongest Correlations 

With our current data set we can address which scaling rela- 
tions are statistical the strongest for nearby galaxies. While 
most of these relations are well known and characterised, we 
discuss several below which are found independent of sam- 
ple selection. These correlations are: T-type with (B — V)o 
colours, T-type with log V max , radius with absolute magni- 
tude, radius with V m ax, (B — V)o with log Vmax, and ab- 
solute magnitude with log V ma x- The relationship between 
these parameters are shown in Figure 10. The best fits be- 
tween these parameters are, 



T = 0.66 ± 0.02 x (B - V) Q 



T = 2.30 ± 0.01 x (log V max J 



0.15±0.02 



-0.06±0.<><>3 



M B = -15.5 ±0.02 x R 



log V max = 1.57 ± 0.001 x R° 



0.1±0.004 



0.16±0.01 



log Vmax = 2.64 ± 0.03 x (B — V)% 



M B = -4.7 ± 0.1 x log V ma x - 9.8 ± 0.2. 

Most of these scalings are already well known, and have 
been described in various contexts previously. Many be- 
come stronger when only considering certain types of galax- 
ies (e.g., ellipticals for the Faber- Jackson, disks for Tully- 
Fisher). We are particularly interested here in the correla- 
tion of T-type with colour and V max (as a proxy for other 
scale features). Figure 10 plots the relationship between 
these properties. As can be seen, later T-types are both 
bluer and have lower internal velocities. These are in par- 
ticlar strong correlations, as shown in Table 3. As we argue 
in this paper, these two correlations are the reason why the 
Hubble sequence has been generally successful at separating 



galaxy types. However, as can be seen there is a large scat- 
ter in these correlations, such that at every T-type there is 
a wide range in colour and V m ax values. This is also true for 
the relationship between (B — V)o and V m ax- This diversity 
cannot be accounted for by errors on classifications, or on 
the measured values. This shows that there are variations in 
colour and scale which is not accounted for by Hubble type. 



6 DISCUSSION 

6.1 Physical Classification 

For the remainder of this paper we address whether the 
results of the first part of this paper can be used to con- 
struct a classification system for galaxies that is physically 
meaningful. First, it is worth asking if we need a physical 
classification system for galaxies. Is it possible that rela- 
tions between observables and (non-morphological) galaxy 
properties are enough to characterise the population? The 
first classification schemes by Wolf (1908), Hubble (1926), 
de Vaucouleurs (1959) and Sandage (1961,1975) are purely 
descriptive, and the idea that a galaxy's classification should 
be based solely on the most obvious features seen in opti- 
cal light persists (Sandage 2005). The reason given for this 
philosophy is that no physical understanding of what drives 
galaxy formation/evolution should be used to create a clas- 
sification system. As Sandage (2005) explains, "A good clas- 
sification can drive the physics, but the physics must not be 
used to drive the classification." 

Sandagc's (2005) argument, which follows Hubble's 
original philosophy, likely resulted from the earliest attempts 
to classify galaxies through the use of theoretical arguments 
by Jeans (1919). Hubble quite rightly avoided theory when 
developing his classification system, yet physical morphology 
does not necessarily rest upon theory, but empirical observa- 
tions. Physical morphology in this sense is not a theoretical 
morphology. If we knew with reasonable certainty what the 
physical effects were that drove the properties of galaxies, a 
classification based on these features would be superior to 
theoretically, or subjectively based systems. 

Several features of galaxies also make them fundamen- 
tally different from other objects which have been success- 
fully classified based solely on their appearance, including 
stars, plants/animals and the chemical elements. Galaxies 
are fundamentally different from these types of objects as 
they are not stable or passive systems, nor do they have 
a well defined life/death as do living organisms and stars. 
Furthermore, elements, plants/animals and most stars spend 
most of their life span in the same form. This is not the case 
for galaxies, where it has been shown that systems of the 
same mass change their appearance in optical light dramat- 
ically in the last 10 Gyr (e.g., van den Bergh et al. 2001; 
Conselice et al. 2005a; Ravindrinath et al. 2006). Galaxies 
are not static systems - they are always changing and evolv- 
ing due to the complex mixtures of stars, gas, dust, dark 
matter, and from interactions with other galaxies. 

Another reason that the physical properties of galaxies 
can be used to devise a classification system is the fact that 
it is not a given that the rest-frame optical light of galaxies 
reveal the most important features of galaxies. For example, 
we now know that the bulk emission from an average galaxy 
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Figure 9. The relation between the projection of galaxies onto the eigenvectors listed in Table 8. The morphological and dynamical 
state of each galaxies is also noted. The red solid points are the ellipticals and SOs, the blue crosses are the disks, and the cyan triangles 
are irregulars. The objects circled and in green are those with an interaction/merger index I > 2. The three projected components, in 
which normal galaxies of different types, as well as interacting/merging galaxies, are well separated, trace the scale (P3), star formation 
(pi), and interacting/merging properties (p2) of galaxies (§5.2). 



is in the mid infrared (e.g., Dale et al. 2005), and if any 
wavelength should be used to classify galaxies it would nat- 
urally be where most of the light is emitted. Fundamentally, 
if we want to classify galaxies through structure, the best 
approach is using a map of the distribution of mass, which 
would be mostly dark matter. A mapping such as this, with 
enough resolution for classification, does not currently ex- 
ist. In blue light, where the Hubble system was and is being 
used, we are examining only the stars in galaxies, and of- 
ten only the brightest ones, which tend to be young. This 
leads to physical classifications, where if we think we under- 
stand the major properties of galaxies, we can utilise these 
for classification. 

Knowing which properties of galaxies are not important 
for classification is just as critical as uncovering the domi- 
nant ones. We can understand why through analogies with 
other sciences. In biological taxonomy, such as separating 
different populations of mammals, features such as hair or 
eye colour, or exact size, are not classification criteria. Sim- 
ilarly, features such as rings, bars and even spiral arms are 
important for the immediate appearance of a galaxy, but 
these are not basic features. Classifications using these cri- 
teria are purely descriptive. Their inclusion as major classi- 
fication criteria are also arbitrary, as other properties, such 
as the steepness of light profiles, the clumpy nature of light 
seen in many galaxies, the spokes and spurs seen in spiral 
arms, etc. were never considered important, or were techni- 
cally infeasible to measure, by early morphologists. 

We can use the results of this paper to construct a 
physically based method for classifying galaxies. Physical 
morphology is not a new idea, and attempts to construct a 
meaningful system for galaxies started with the work of Mor- 
gan (1958; 1959) who attempted to correlate "the forms of 
certain galaxies and their stellar content as estimated from 
composite spectra" from Morgan & Mayall (1957). Later at- 
tempts include van den Bergh's (1960) study on the corre- 



lation between spiral arm structure and intrinsic brightness. 
Modern studies have attempted to classify ellipticals by their 
structures (Kormendy & Bender 1996), and at using inter- 
acting and star formation properties to classify all galaxies 
(e.g., Conselice 1997; Conselice et al. 2000a; Conselice et 
al. 2000b; Bershady et al. 2000). As we argued in §5, star 
formation, interactions/mergers, and the masses of galaxies 
are the three most fundamental parameters that can be used 
together to form a new physical galaxy classification. 

6.1.1 Mass 

We concluded in §5 that the underlying mass (scale) of a 
galaxy is a fundamental property, which is responsible for 
the major component in our PCA analyses. While few would 
argue that the mass of a galaxy is not a fundamental cri- 
terion for classification, it is not obvious that it is possible 
to classify galaxies based on their mass, without direct mea- 
surements. 

However, there are structural indicators in optical light 
which can be used as a proxy for the underlying mass in a 
galaxy. As shown by e.g., Graham et al. (2001) and Con- 
selice (2003) the concentration of a galaxy's light, measured 
through either parametrised fits to light distributions, or 
from direct measurements using a concentration index (Ber- 
shady et al. 2000), broadly correlate with the stellar and 
kinematic masses of galaxies. This is true at low redshift, as 
well as out to z ~ 1 (Conselice et al. 2005a). 

The degree of apparent concentration of stellar light in 
a galaxy has been used (even if not explictly stated) since 
Hubble (1926) as the main axis for nearly all galaxy classi- 
fication systems. As shown in e.g., Conselice (2003) galaxy 
concentration relates to the bulge to disk ratio and Hub- 
ble type, although this breaks down for dwarfs. Further- 
more, light concentration has continued as the main axis in 
all galaxy classification systems, such as in de Vaucouleurs 
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Figure 10. The strongest correlations found between the properties of the 22,000+ galaxies studied in this paper. Most of these 
relationships are already well known, and described for morphological subsets, such as the Tully-Fisher relation, and the size-magnitude 
relation. For our purposes, the most interesting trends arc those of T-type with colour and V max , where considerable scatter can be seen. 
This is also true for the correlation between colour and V ma x- The scatter in these correlations is larger than the observational errors, 
and suggests that while morphology, colour, and scale all correlate in nearby galaxies, they are still independent features. 
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Figure 11. The three dimensional classification volume proposed in this paper to account for all known galaxies. The three parameters 
used are the mass, recent star formation and interaction properties of galaxies. Several well known galaxy types are labelled in their 
appropriate locations in this space. At the top of the classification cube are found galaxies with no recent star formation. Across the 
mass axis we progress from dEs to SOs to giant ellipticals. Any galaxies with a recent interaction/merger would be found further into the 
third dimension. As we go further into the cube within the y-axis, we find systems with a higher contribution of recent star formation to 
the stellar mass. Low-surface brightness disks would be found in a region of this space where the star formation is modestly high, and 
interaction history low. Merger induced starbursts are located in the high star formation, high interaction part of the diagram, while 
irregulars and dls are located towards the bottom right of the cube, as these are typically lower mass star forming galaxies. The Hubble 
sequence is a one dimension projection of the star formation and interaction axes onto the mass sequence. This classification necessarily 
replaces idealised Platonic Hubble types by idealised parameters that are not always easy to retrieve, but in principle are measurable for 
all observable galaxies. 



(1959) and Morgan (1959), and continues as a quantitative 
system through fitting bulge to disk ratios. 



While bright Hubble types generally correlate with the 
masses of galaxies, they do not account for other galaxy 
properties needed for a full classification. The fact that late- 
type galaxies, which generally have low: masses, bulge to 
disk ratios, and light concentrations, are blue is a result of 
evolutionary processes, rather than an inherent success of 
the classification. It is a coincidence that Hubble types cor- 
relate well with colour in the nearby universe, as this is not 
the case at higher redshifts where morphologically classified 
ellipticals are often blue in colour and star forming (e.g., 
Stanford et al. 2004; Bundy et al. 2005). This correlation 
breaks down even more at z > 1, where the Hubble sequence 
itself begins to dissipate. 



6.1.2 Star Formation 

As shown in §5 we find a morphological dichotomy between 
galaxies that are dominated by recent star formation and 
those that are not. Our principal component analyses re- 
veals that colour is the major contribution to a significant 
principal component, which is not accounted for solely by 
the Hubble type, as galaxies at a given morphology type 
have a range of colours (§5.1, Figure 10). Therefore, if we 
knew the current star formation rate of a galaxy it would 
be a very powerful classification parameter. Typically the 
current star formation rate is derived from either ultravi- 
olet, Ha or far- infrared fluxes (Kennicutt 1998), although 
these measurements are typically time-consuming and hard 
to obtain for large samples. It is therefore desirable to use a 
method that can estimate the star formation properties of 
galaxies based on optical imaging. For most galaxies, this 
can be done by combining colour information along with 
high-frequency structural features produced by young stars 



18 Christopher J. Conselice 



(Isserstedt & Schindler 1986; Conselice et al. 2000a; Con- 
selice 2003). The dumpiness index, S, (Conselice 2003) is 
a good measure of the unobscurbed current star formation 
rate in galaxies, and can therefore be used as an indicator. 

6.1.3 Galaxy Mergers/Interactions 

Like the scale of a galaxy and its current star formation rate, 
galaxy mergers and interactions are an important aspect for 
driving galaxy evolution. Although mergers are rarely found 
in the local galaxy population, they have strong effects on 
the evolution of galaxies and their properties. Ellipticals and 
high surface brightness bulges are likely the result of major 
mergers, and other accretion processes. In reality, all merger 
events for any one galaxy are nearly impossible to trace with 
certainty, especially old mergers that have gone through vi- 
olent relaxation, after which the original orbital information 
is lost. Recent interactions between galaxies however will al- 
most always have obvious signatures that will change the 
way galaxies appear, as well as their internal evolution, as 
has been shown in various methodologies. 
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Figure 12. The system in Figure 11 projected onto the star for- 
mation/mass plane. The Hubble sequence is labelled, and includes 
galaxies with a range of masses, modest star formation with no 
distinction made for interactions. Ellipticals, SOs, dEs and dwarf 
spheroidals are located at the bottom of the star formation se- 
quence. Massive star bursts are those galaxies located at the high 
star formation area of this projection. 



6.2 A New Classification System: The 
Classification Cube 

Figure 11 shows how galaxies of various types fit into a 3- 
dimensional parameter space defined by galaxy mass, recent 
star formation, and interactions. On this diagram, the x-axis 
signifies the mass, while the y-axis represents the recent star 
formation, and the z-axis represents the degree of recent 
interactions with other galaxies. 

As opposed to previous classification systems, the one 
presented in Figure 11 can account for all the major galaxy 
types described in §2. Classical ellipticals, with little sign of 
recent interactions or star formation will occupy the high- 
mass, low star formation and low interaction corner. The top 
of this classification volume also contains the SOs, dwarf el- 
lipticals/spheroidals. SOs have little evidence for recent star 
formation, but are lower-mass systems than ellipticals, and 
they occupy the same star formation and interaction space 
as the ellipticals. dEs are systems with low-mass, and as 
such they occupy the far end of the mass sequence for these 
galaxies. The dwarf irregulars have a similar internal struc- 
ture and mass to the dEs, except they contain star forma- 
tion, and thus fall in the same mass and interaction space, 
but are separated in star formation space. 

This system is not completely new and various aspects 
of it have been published in different forms. For example, van 
den Bergh (1976) proposed a revision of the Hubble system 
where SO galaxies occupy a parallel sequence to the spirals, 
instead of lodged between the ellipticals and Sa galaxies as in 
the original Hubble sequence, van den Bergh's (1976) classi- 
fication also contains an anaemic galaxy sequence for objects 
with little gas. This system also classifies galaxies in terms 
of their bulge to disk ratios, making it very similar in some 
ways to the star formation/mass plane presented here (see 
Figure 12). 

Low surface brightness galaxies, as discussed earlier, 
have little evidence for companions or signs of interactions 
with other galaxies. They also have a very gradual star for- 
mation history that is still ongoing. As such they occupy the 
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Figure 13. Projection of the star formation/interaction plane 
of Figure 11. Several classes of galaxies separate along this plane 
perpendicular to the mass sequence. Low surface brightness galax- 
ies (LSBs) are those with modest star formation and very little 
recent interaction history. Ellipticals and SOs are objects with 
relatively little star formation or interactions, and no recent in- 
teractions/mergers. Several other galaxy types, such as distorted 
cluster galaxies undergoing little star formation due to gas re- 
moval by the cluster, dry-mergers, and starbursts induced by ma- 
jor mergers, are also labelled. 



same mass sequence as the spirals, but are located on the 
low interaction side of the interaction sequence. 

The dichotomy between ellipticals that appear boxy and 
disky can also be represented by this system (Kormendy & 
Bender 1996), since these galaxies naturally fall along differ- 
ent regions of interaction and mass space. Disky ellipticals, 
with rotation, would have a lower mass, while the degree of 
boxiness correlates with the interaction history. 

Other galaxies can also be classified in this system. Star- 
bursts induced by major mergers, are placed at the high 
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end of the interaction and star formation sequence. Galaxies 
undergoing violent relaxation, such as Arp 220, containing 
spheroidal components and step light profiles (Wright et al. 
1990), would be classified towards the high-mass end of the 
mass sequence. 

This system is constructed based on the results pre- 
sented in the earlier part of this paper. Measuring the star 
formation, merging, and masses (either stellar or dynami- 
cal) properties of a galaxy are not trivial, but are likely to 
become easier in the future. We are therefore not setting out 
in this section a stringent method for classifying galaxies by 
these properties. One could do these classifications directly 
by measuring the star formation rate through emission lines, 
masses through dynamics or spectral energy distributions, 
and mergers/interactions through an asymmetry index or a 
dynamical indicator. We discuss in §6.3 a possible way to 
use the apparent morphologies of galaxies in optical light to 
carry out these classifications. Future work will reveal better 
ways of carrying out these measurements, which can then be 
applied for classification according to Figure 11. 

6.2.1 Other Systems as Projections 

Other morphological systems, including physical morpholog- 
ical ones, can be viewed as projections of this system. The 
Hubble sequence, to first order, is a one dimension projec- 
tion with the star formation and interactions collapsed onto 
the mass sequence. As discussed before, the Hubble sequence 
does not take into account these features, thus there exists 
a wide range of properties at each Hubble type (Roberts & 
Haynes 1994; Kennicutt 1998). 

The colour-asymmetry diagram is a projection, in part, 
of the star formation, interaction sequence (Figure 13). This 
is not an exact projection since the asymmetry parameter 
measures interactions, as well as star formation (Conselice 
et al. 2000a). To first order the (B — V)o colour measures 
the star formation history of a galaxy. Galaxies with purely 
old stars have very red colours, while those with recent star 
formation are bluer. The concentration-asymmetry diagram 
used by Bershady et al. (2000) and others, is furthermore a 
projection of the mass-interaction sequence. More recently, 
Driver et al. (2006) argue that light concentration and colour 
are fundamental properties that separate galaxies into early 
and late types, which is again a two-dimensional projection 
of this new system. Likewise, Blanton et al. (2005) sug- 
gest that colour and and luminosity are the most important 
galaxy properties, which is also shown in our analysis. 

6.2.2 Idealised Classifications 

When Hubble invented his classification, he did so with spe- 
cific examples of galaxies, however Hubble types are now 
thought to characterise ideal types. There are few galaxies 
that fit this ideal, nor is any one galaxy exactly similar to 
any other. This is due to the chaotic interaction and star 
formation history of galaxies that are, at present, not easily 
retrievable. 

The classification system presented here replaces ideal 
galaxies with parameters that are in some ways idealised 
measurements. It is not trivial to determine the star for- 
mation history of a galaxy, or even accurately measure its 



current star formation rate. Star formation does have a mor- 
phological appearance in most galaxies, the exception being 
systems heavily enshrouded by dust. The system presented 
here can still account for those galaxies, as it is very general. 
The method for classifying the star formation properties of 
a galaxy does not have to be based on morphology, but can 
also be done through other properties, such as fluxes at var- 
ious wavelengths. 

The same is true for measuring the merger history of 
galaxies. While many galaxies have some clues left over from 
merger events, these clues can be quite diverse, requiring 
various methods to retrieve them. Structural methods, such 
as the asymmetry CAS parameter (Conselice et al. 2000a,b), 
are useful for determining the presence of, and quantifying, 
current major mergers. Determining the history of past in- 
teractions is more difficult. One way is to investigate the 
outer parts of galaxies where the dynamical time scales are 
long, and where apparently normal galaxies show evidence 
for past merger activity (Schweizer & Seitzer 1988). Other 
methods include searching for distorted HI halos surround- 
ing galaxies, galaxy core depletion from merging black holes 
(Graham 2004), and finding distorted rotation curves (Swa- 
ters et al. 1999) or line profiles (Haynes et al. 1998; Conselice 
et al. 2000b). These features are potentially all remnants of 
accretion, merging or interactions with other galaxies, and 
can in principle be used to quantify the interaction history 
of galaxies. 

6.2.3 Unidentified Galaxies? 

The system presented here can classify all known galaxies, 
but it also has classification 'corners' where no nearby galax- 
ies are known to exist, or have not been identified. One type 
are massive, low interaction high star formation systems. 
There are plenty of galaxies that have these features at high 
redshift. These are likely ellipticals or disks in formation that 
will evolve into spheroids. For example, systems with high 
mass, high star formation, and little merging could origi- 
nate from gas accreted onto elliptical-like systems. Other 
possible galaxies not yet identified are those with various 
mass ranges, but low star formation and high interaction 
strength. Perhaps the population of distorted, but non-star 
forming galaxies, or compact dwarfs found in clusters (Con- 
selice & Gallagher 1999; Drinkwater et al. 2003) are exam- 
ples of these, where the gas has been removed by the in- 
tracluster medium. Another example are the 'dry' mergers 
between spheroids with little gas (e.g., Hernandez- Toledo et 
al. 2006) . Other possible galaxy populations can be hypoth- 
esised by studying Figure 11. 

6.3 An Optical View of Galaxy Classification 

We have argued that the three parameters: mass, star for- 
mation and interactions can be used together to classify and 
trace the evolution of galaxies. These parameters have a di- 
rect correlation with the semi-irreversible processes of mass 
assembly or removal, and star formation where gas mass 
is converted into stellar mass by some conversion efficiency 
q. Therefore, in principle these features can be measured 
through stellar light distributions, which has been the pri- 
mary method for classifying galaxies for nearly a century. 
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We show in this section that this is possible. The follow- 
ing description follows closely the more detailed calculation 
given in appendix B of Conselice (2003) . 

We can write the initial baryonic mass of a galaxy after 
its initial formation, (Mq,b) as, 



M ,b = M ,s +M ,g, 



(1) 



where Mo,s and Mo,g are the initial stellar and gas masses. 
Over time, mass will accrete or be removed from a galaxy 
due to interactions and mergers. At a time T, the amount of 
mass added to or removed from the system, M'(T) is given 
by, 



M'(T)= [ M(t)dt= ( M s (t)dt + f M G (t)dt, 
Jo Jo Jo 



(2) 



with M the mass accretion, or removal, rate as a function 
of time, and T is the age of the galaxy when observed. We 
can further divide the mass flux into stellar Ms and gas 
Mg components. Some gas will be converted into stars at 
each time, t, by a conversion rate, a(t), such that the mass 
converted into stars, Ms, is given by 
r-T 



M S (T)= / a{t)M G (t)dt, (3) 
Jo 

where a(t) and Mg are the gas conversion rate, and gas mass 
as a function of time. We can then write the total stellar and 
gas masses as a function of time by: 

M s (t) = M ,s+ / (M s (t) + a(t)M G (t))dt, (4) 
Jo 

M G (i) =M ,g+/ (M G (t) - a(t)M G (t)) dt. (5) 
Jo 

The three axes of our system (Figure 11) can be represented 
by three terms in these equations. An interaction parameter 
can be represented by Mg (t) , and a star formation parameter 
by Ms. We can also create dimensionless parameters that in 
principal correlate with measurable stellar light features. We 
can write the integrated interaction, I(t) and star formation 
F(t) strengths weighted by the final mass, and time as: 



lit) 



F{t) 



t \M s (t)\ 
T M S (T) 



dt 



t a(t)M G (t) 
T M S (T) 



dt. 



(6) 



(7) 



These two ideal morphological parameters are quantities 
that should correlate with the appearance of a galaxy, and 
are further elaborated in appendix B of Conselice (2003). 
They are weighted by time since older star formation and 
interactions will have less of an effect on the appearance of 
a galaxy. Likewise, more massive galaxies will have appear- 
ances that, in general, will be less affected by relatively lower 
amounts of gas converted into stars, and stars lost or gained 
from galaxy interactions. 

What this section shows is that the main galaxy forma- 
tion processes - star formation and interactions/mergers can 
be traced through optical light. Typically, a galaxy which 
has undergone recent star formation will contain clumpy 
light, occupied by clusters of young stars. Likewise, a galaxy 
undergoing a merger or interaction will appear lopsided 



and/or distorted in some manner. The mass of a galaxy 
often correlates with the light concentration, which can be 
estimated by eye. Therefore, by estimating the degree of: 
disturbances, clumpy structures, and how concentrated a 
galaxy is (i.e., bulge vs. disk structures) it is possible to 
classify galaxies into this new system using optical light. A 
detailed exploration of this idea in terms of eye-ball classi- 
fications is however beyond the scope of this current paper. 
Quantitatively, this approach has been successfully used to 
classify both nearby (Conselice et al. 2002; Conselice 2003) 
and high redshift galaxy populations (Conselice et al. 2003b; 
Conselice et al. 2004; Grogin et al. 2005; Papovich et al. 
2005; Lehmer et al. 2005; Conselice et al. 2005c; Lotz et al. 
2006) within the CAS system. 



7 SUMMARY AND CONCLUSIONS 

This papers attempts to answer a few very basic questions 
about galaxies and their properties. The first is how galax- 
ies of different morphological types are related to each other, 
including the relationship between various major evolution- 
ary processes. Another question we address is whether or 
not the structural appearance of a galaxy contains enough 
information to offer this classification. To answer the first 
question, which constitutes the bulk of the paper, we car- 
ried out a series of statistical analyses, including principal 
component analyses, of fourteen properties in 22,121 galax- 
ies. The major results of our investigations are: 

I. Spearman correlation analysis show that most measur- 
able properties of galaxies are correlated with each other, 
and that Hubble types correlate strongest with stellar mass 
and (U — B)o and (B — V)o colour. Scale features (M*, 
Mb, radius, HI gas content, far-infrared luminosity) corre- 
late strongly with each other independent of morphological, 
colour, or luminosity selection. Furthermore, the strong cor- 
relation between M* and colour is independent of Hubble 
type. 

II. We carry out several principal component analyses 
(PCA) to determine how many PCA eigenvectors are needed 
to describe the major properties of all nearby galaxies. Our 
conclusion from this is that the major eigenvector is domi- 
nated by the scale, or mass, of a galaxy, while the next two 
major eigenvectors are dominated by the star formation, and 
galaxy interactions/mergers, respectively. 

III. Based on points I and II, we conclude that the three 
major processes that should be used to classify a galaxy 
are: the scale or mass, interactions and mergers with other 
galaxies, and recent star-formation. 

IV. We argue that these three features (mass, star forma- 
tion and interactions) are revealed through the stellar mass 
distributions within galaxies. Thus the optical structure of a 
galaxy can be used, in principle, to accurately classify galax- 
ies into a physically meaningful system. One way this can be 
done is through the CAS parameters described in Conselice 
(2003). 

V. Based on our analysis we present a morphological break 
down of z ~ galaxies, including the fraction and num- 
ber densities of different Hubble types (ellipticals, spirals, 
irregulars) at different luminosities. We find that the bulk 
of nearby bright galaxies with Mb < —20 are spirals (69%), 
while 19% are elliptical, 10% SO and 2% irregular. We also 
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investigate which types of galaxies are barred, and the num- 
ber densities of barred galaxies in the nearby universe. These 
are intended to be useful for higher redshift comparisons, 
which are just now becoming possible. 

In summary, this paper presents a new basic methodol- 
ogy for understanding and classifying galaxies. It is a phys- 
ical system that uses the quantifiable parameters of mass, 
star formation, and interactions, and is not based on galaxy 
morphology. The system can account for all known galaxy 
types, and is a naturally outgrowth of detailed studies of 
galaxies over the past century. Future work on this classifi- 
cation scheme will include understanding the origin of the 
mass sequence, and why it correlates at z ~ with star 
formation and bulge to disk ratios. Furthermore, our new 
classification system is by no means finalised. In the future, 
we will likely be able to measure galaxy masses, star forma- 
tion rates, and ongoing interactions with a high accuracy. 
Whether these measurements are done morphologically or 
not is irrelevant, as these three properties measured in any 
way can be used to classify galaxies. 
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